Oxide superlattices with alternating p and n interfaces 
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The physics of oxide superlattices is considered for pristine (001) multilayers of the band insulators 
LaAlOa and SrTiOa with alternating p and n interfaces. A model of charged capacitor plates offers 
a simple paradigm to understand their dielectric properties and the insulator to metal transition 
(IMT) at interfaces with increasing layer thickness. The model is supported by first-principles results 
based on density-functional theory. The charge at insulating interfaces is argued and found to be 
as predicted from the formal ionic charges, not populations. Different relative layer thicknesses 
produce a spontaneous polarization of the system, and allow manipulation of the interfacial electron 
gas. Large piezoresistance effects can be obtained from the sensitivity of the IMT to lateral strain. 
Carrier densities are found to be ideal for exciton condensation. 

PACS numbers: 73.21. Cd, 71.30.-|-h, 73.20.-r 



I. INTRODUCTION 

Recent technological advances have enabled the fab- 
rication of high quality oxide multilayers, revealing a 
wealth of fascinating new physics. One striking exam- 
ple is the LaAlOs/SrTiOa system (LAO/STO). In 2004, 
Ohtomo and Hwan^ discovered that the interface be- 
tween these two perovskite band insulators can be con- 
ducting, depending on the termination of both mate- 
rials. Many experimental and theore tical studies have 
followed^™ g l ^ l '^ l ^ | i" l " | i" | i^ | i^ | i^ | i^ l " l. Charge compen- 
sation at the interface of thick layers is required to avoid 
the so-called polar catastrophe. It arises from the di- 
verging electrostatic potential caused by the net electric 
charge at these interfaces resulting from the fact that the 
(001) monolayers of STO are neutral while the ones of 
LAO are charged. The Ti02-LaO interface (termed n) 
needs 0.5 electrons per two-dimensional (2D) unit cell, 
and the SrO-A102 interface (p) 0.5 holes to avoid this 
polar catastrophe. Such numbers are based on formal 
ionic charges, i.e., Ti+'*, Sr+^, La+^, Al+^, and O^^, al- 
though it is well known that the charge distribution in 
these materials is far from being so ideally ionic. For- 
mal charges are often downscaled by so-called covalency 
parameters that aim to describe more realistic charges^. 

These compensating electrons and holes are confined 
to the interface regions, bu t hig hly mobile in 2D, defining 
2D electron and hole gased^l^ Huijben et alW (se^^^ for 
the theory) studied a system with one p and one n in- 
terface, finding 2D conduction beyond a critical interface 
separation of five unit cells. Characteristics like carrier 
mobility depend on the carrier density, which grows with 
separation beyond the critical valuelSj. Multilayers with 
both p and n interfaces with such a control of carrier den- 
sity provide promising systems for obtaining the so-far 
elusive excitonic insulatoil^. Recent experiments have 
demonstrated that oxygen vacancies can be controlled, 
while sti ll all owing metallic conduction attributed to a 
2D layer ISESl. 



This paper focuses on superlattices with alternating 
p and n interfaces. We consider pristine syste ms, wi th 
no point defects (the effect of oxygen vacancie^^^ESl jg 
discussed at the end). Considering an equal thickness 
for both materials, an IMT is found by our calculations 
based on density-functional theory (DFT) for a thickness 
of just over eight unit cells. The electrostatic potential 
obtained in the calculations agrees remarkably well with 
a simple model of capacitor plates, giving an almost con- 
stant field of opposite sign in both materials, which does 
not change with increasing thickness until the potential 
drop coincides with the band gap. At that critical thick- 
ness electrons transfer from the p to the n interface, mak- 
ing them conducting, and pinning the potential drop. 



II. METHOD 

The DFT calculations were done using t he loc al den- 
sity approximation^^ and the Siesta methocP^Ml. Norm- 
conserving pseudopotential^^l were used, considering 
normal cores for O and Al, while semi-core electrons 
were included in the valence for La(5s5p), Sr(4s4p) and 
Ti(3s3p). Double-^ polarized bases were used for valence 
electrons'^. Integrals in real space were performed on a 
mesh of 250 Ry cutofl'^^, while Brillouin zone integrations 
were done on a fc-mesh of 30 A cutofpSl. Four unit cells 
of STO in its ideal perovskite structure were layered on 
top of four of LAO, the 4/4 superlattice containing both 
Ti02-LaO and SrO-A102 terminations. The lateral lat- 
tice parameter was set to the theoretical average for the 
two materials. The cell size was relaxed perpendicular 
to the interface (z) along with the atomic positions un- 
til the forces were below 15 meV/A. Samples for 8/8 and 
12 /12 unit cells were equal ly prepared. Other DFT stud- 
ies [MIlITaiTITlIilIiliI^ have focused on single inter- 
faces of either kind, on arrays of repeated in terface s, or 
on alternating p and n interfaces. The lattei'^'^'^ have 
considered four unit cells or less, too thin to observe the 
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FIG. 1: Band structures of three LaAlOs/SrTiOs superlat- 
tices with (a) 4/4, (b) 8/8, and (c) 12/12 unit-cell thicknesses. 



physics discussed here. 



electrons are confined to the Ti02 plane at the n in- 
terface, and around T. The 2D electron gas (2DEG) is 
dominated by a Ti 3d character. It is nicely parabolic 
with an effective mass of 0.4 rrie- The several parallel 
sub-bands at this interface correspond to excitations un- 
der the effective confining potential for electrons along z 
(Fig. 4ap. 

Quantitative predictions for the IMT are biased by 
the known band-gap problem of Kohn-Sham fermiona^Sl. 
Constrained DFT calculations as used for charge-transfer 
systemd^can be used for ours, but are beyond the scope 
of this work, its focus being the elucidation of the main 
mechanisms at play. In addition, it has been argued that, 
being LaTiOa a system with highly correlated electrons 
in the Ti Sd band, methods addressing strong correla- 
tions are needed for our system. Note that this is not 
the case, however, since our IMT is in the limit of zero 
occupation of the 3c? band, hardly a correlated system, 
unlike LaTiOa, which has one Sd electron per Ti atom. 

The deviations from the zig-zag behavior of the poten- 
tial are small in the scale of its amplitude, allowing the 
definition of a net electric field through each material. Its 
magnitude is nearly equal in both materials despite their 
different dielectric response. Such response is illustrated 
in Figure 2 by the z-splitting of cations and anions at 
each layer, which is much larger for STO as expected. 
The fields obtained are 57.3 mV/A and 57.1 mV/A for 
the 4/4 and 8/8 superlattices, while the 12/12 sees it 
reduced to 37.8 mV/A due to the partial charge back- 
transfer. The field values correlate with the difference in 
sub-band separation seen in Figure 1 . An estimate of the 
2DEG width W is obtained from W ~ Ezpe/£, i-e. the 
zero-point energy (ZPE) for the confining potential, over 
its slope, the electric field. Taking for 8/8, £ = 57 mV/A 
and EzpE ~ 0.2 eV gives ~ 4 A. 



III. RESULTS AND DISCUSSION 
A. Superlattice calculations 

The band structures of the three superlattices are pre- 
sented in Figure 1. The band gap is indirect and reduces 
with interface separation, closing for the 12/12 system, 
with holes and electrons separated in reciprocal space. 
The behavior in real space is shown in Figure 2 for the 
12/12 case using the density of states projected onto each 
bilayer. The electrostatic potential is plotted alongside. 

The physics of the problem is apparent in this figure. 
The net electric charge of chemical origin at the interfaces 
establish electric fields among them, defining a zig-zag 
potential, which is closely followed by the band structure 
of both materials. The IMT occurs when the amplitude 
of the zig-zag is larger than the band gap (the valence 
band offsets are small in this scale, 0.2 eV for the p and 
0.0 eV for the n interfaces). At this point holes appear 
confined in z around the AIO2 plane of the p interface, 
and around the M point in the 2D Brillouin zone, while 



B. Model of charged plates 

The physics described can be further analyzed with a 
simple model. The insulating system is modeled by a se- 
quence of capacitor plates, one per interface, separated 
by dielectric material, as sketched in Figure [3j Each plate 
has a planar charge density of chemical origin Uc (the 0.5e 
per interface unit cell described above), positive at the 
n interface, negative for the p. Using Gauss's law, peri- 
odic boundary conditions (PBC), and equal thicknesses, 
there is a uniform electric field of magnitude £^ pointing 
outwards from the n interface that satisfies 

- Plao - ^STO = 2eof °, (1) 

with Plao E^nd PsTO the magnitude of the respective 
polarizations of both materials under the field, and eg the 
dielectric permittivity of vacuum. The left hand side is 
what is indicated by a net in Figure [3] The behavior at the 
p interface is exactly opposite, with the field now pointing 
towards the plate. The fields have equal magnitude in 
both materials even if both polarizations are different. 
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FIG. 2; Centre: Density of states projected on unit-cell bilayers for the 12/12 superlattice. Left: AveragecP^ electrostatic 
potential energy for electrons as a function of z. Right: The anion-cation splitting of the Ti02 planes in SrTiOs and LaO 
planes in LaAlOa (larger splitting in each material), with 5z = Zcation — ^anion- 



Finding £° requires knowing P{£) beyond linear re- 
sponse, at least for STO. Under the strain conditions im- 
posed by our geometry, the bulk o f STO presents a spon- 
taneous polarization along JSoEl]^ which allows a simpler 
modeling than for unstrained STO, and offers an upper 
bound in the response, which corresponds to an upper 
bound to the critical thickness for the IMT. We take 
^LAO = eoXLAcS and Psto = -Psto + '^oXsto^- Defining 
K = 2 + XSTO + XlaOj equation 1 becomes 

= (a, - PsTo)/eo^. (2) 

We obtain Pgxo = 0.309 C/m^ using the Berry phase 
approach'^^ for bulk STO with the same strain condi- 
tions as in the multilayer. The lattice contribution to 
both susceptibilities is computed as in Ref. from the 
phonons and the Born effective charges, which are ob- 
tained by finite differences of the forces and the polariza- 
tion, respectiveljESl. The computed susceptibilities are 
xf3(LA0)= 12.2 and xf3(ST0)= 24.7. To these we 
add the electronic contribution taking it from e°°(STO)= 
S.l^and e°°{LAO)=4.7i^. Using these bulk quantities 
and (Tc = 0.5e per interface unit cell, an electric field of 
57.4 mV/A is obtained, in excellent agreement with the 
superlattice DFT results. 



C. Chemical charge at interfaces 

This agreement indicates that the physically meaning- 
ful charge CTc is the one predicted by formal ionic charges, 
rather than populations. The reason is the same as 
for dopants in semiconductors, where irrespective of the 
charge distribution around the dopant, the net charge is 
the result of counting core charges on one hand and elec- 
trons in the valence band on the other. Take the example 
of phosphorous as a donor in bulk silicon. Its -f 5 core is 
surrounded by the same set of bonds that surrounded 
the Si -1-4 core it substitutes. The electron-pair clouds of 
these bonds are deformed (polarized), of course, but the 
number of electrons remains. The net effect is that the 
P donor creates an attractive potential for electrons that 
corresponds to a charge of +e, not an effective charge 
of any kind, the deformation of the bonds around being 
described by the dielectric constant of bulk silicon. 

These arguments were generalized to interfaces three 
decades agcpSI. The valence-band electron counting was 
performed in terms of bands, but it can also be done 
using Lewis pairs or localized Wannier states. In LAO 
and STO the valence band corresponds to four Wannier 
states per oxygen, which gives exactly the same numbers 
as when using formal charges. The local charge neutrality 
picture^^ can then be used as the easiest way to see the 



FIG. 3: (Color online) Capacitor-plates model of the 
LAO/STO superlattice. Above, the plates are indicated by 
the thinner bands, and ctc indicates the charge of chemi- 
cal origin attached to each, which is equal but of opposite 
sign for alternating interfaces. The box around the central 
plate indicates the surfaces for integration of Gauss's law. 
The lower panel shows the net electrostatic potential V for 
the system, which can be seen as arising from plates with 
Tnet = cTc — Plao — -PsTO, the latter Plao and Psto referring 
to the magnitude of the respective polarizations (notice that 
V has opposite sign to what shown in the left panel of Fig. 2, 
which is qV). 

interface charging, but using formal charges, not effective 
ones. Whether atoms are chemically more or less ionic 
depends on the shape and displacement of the Wanniers 
from the O atom towards the cations, but is irrelevant 
here. What is relevant is the number of electron pairs 
and the fact that they localize over lengths much smaller 
than the interface separation. 

D. Insulator to metal transition 

The IMT occurs at the critical thickness of dc = A/f °, 
where A is the gap (STO's), including the ZPE for the 
confining potential at both interfaces (Fig. 4a). Tak- 
ing the calculated £'^, our DFT gap for strained STO 
(1.78eV), and both ZPEs as - 0.2 eV, a 4 = 38 A is ob- 
tained, ~ 10 unit cells. This value is sensitive to lateral 
strain, since P{£) depends on strain (especially for STO, 
so close to a ferroelectric instability). Indeed, repeating 
the DFT calculation for 8/8 under the lateral lattice pa- 
rameter of bulk STO gives conducting interfaces. This 
effect can be used to sense applied strain on a sample 
tuned to be close to the IMT (a piezoresistance device). 

Beyond the IMT, a constant density of states for holes 
and electrons is assumed, modeling 2D parabolic bands. 
Considering electron transfer from the p to the n inter- 
face, equation (1) becomes 

- Plao - ^'sto - {Sd - A)D = 2eoS (3) 

where d is the interface separation, £ the modified field, 
and D the density of states (taking equal effective masses 



FIG. 4: (Color online) Model band gap versus z for a su- 
perlattice with equal thicknesses (a), the thin film sample of 
reiW (b), and a superlattice with different thicknesses with 
open (c) and periodic (d) boundary conditions. In (a) and 
(c) the arrow indicates the charge transfer leading to electron 
(red) and hole (blue) 2D gases. In both cases the thicknesses 
depicted give the onset of the charge transfer, corresponding 
to the insulator to metal transition. 

for electrons and holes). Proceeding as for equation 2, 
£^ia, + DA- PsTo)/(eo'« + Dd). (4) 

The electric field vanishes as 1/d for large separations, the 
charge transferred tending to compensate the chemical 
charge. In our case, DA ^ (Tc, and thus £/£^ ^ dc/d, 
i.e. the voltage drop is essentially pinned by the gap. 

E. Extending the model to other systems 

Consider now the model for other systems. Take a 
superlattice with thicker layers for LAO than for STO 
(Fig. 4c). It is analogous to a frozen (111) LO phonon 
in NaCl: the different separations of the charged planes 
gives a net electric field. PBC remove it as if putting the 
system between shorted metal plates (Fig. 4d). It is as a 
ferroelectric, except that switching the spontaneous po- 
larization demands changing the thicknesses. Note that 
with PBC and thinner STO, fsTO > ^lao- Without 
shorted plates the system is unstable to the appearance 
of 2DEGs (Fig. 4c). 

The dc found experimentally on the non-repeated sys- 
tem of ref.l^is much smaller than our superlattice result. 
They have one LAO layer interfacing n to the STO sub- 
strate and p to an STO overlayer (Fig. 4b). The field is 
zero in the substrate and overlayer, the p and n inter- 
faces defining a capacitor only screened by LAO. Thus, 
£ = CTc/eo(l + Xlao), more than twice as large as £^ 
of equation 2, giving the observed smaller dc- A similar 
model has been very recently described for one single n 
interfac^^. 

IV. FINAL REMARKS AND SUMMARY 

The electron-hole interactions among 2D gases should 
establish excitons. Achievable exciton densities are fa- 
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vorable for exciton condensation, the optimal density be- 
in^nopt = (ttoq)"^, with Oq « {aBe/rn'Y+cP (Bohr ra- 
dius, dielectric constant and effective mass, respectively). 
For our system Uopt ~ 0.002e/cell, well within range 
(ni2/i2 = 0.073e/cell). 

The IMT is affected by the presence of O vacancies. 
Each donates two electrons to the n interface giving rise 
to a 2DEG. Their appearance and location depends on 
their stability and on kinetic effects like electromigration 
and sample-growth history. Taking stability arguments 
only, the IMT via vacancies would happen when S'^d > 
Ho/2, with £^ as in equation 2 and /io the formation 
energy of an O vacancy at the p interface, which depends 
on the oxygen chemical potential at growth conditions. 



In summary, the agreement between DFT and the 
model shows that it contains the relevant physics of these 
superlattices, pointing to new science and applications by 
changing relative thicknesses, substrates, and stress. 
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